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I 1 Introduction 

Oh! 
(D . 
^ , In this talk, we will discuss high-order corrections to heavy quark production 

in QCD. To get a better confidence in our ability to predict cross sections 

k>«( ; for various heavy quark reactions, including open fiavor and quarkonia pro- 



a 



Abstract 

We review some techniques of resummation apphed to heavy quark 



^j ■ duction, we need to gain control over a class of corrections associated with 

what is sometimes called partonic threshold. This talk will describe the na- 
ture of these corrections, and report on the progress of the past few years 
[H, |, ^ |, U, H in their resummation to all orders in perturbation theory. 

Corrections due to resummation may turn out to be small, in which case 
our confidence in low-order perturbative cross sections should increase, or, 
they may turn out to be large, and may afford sensitivity to QCD dynamics 
in a regime where all orders of perturbation theory are relevant. It is possi- 
ble that both scenarios apply in different cross sections, and/or in different 
kinematic ranges of the same cross section. 

^Talk presented at Workshop on Physics With Electron Polarized Ion Collider, EPIC 
'99, April 8-11, 1999, Indiana University Cyclotron Facility. 



2 Threshold Resummation for Hard Inclu- 
sive Scattering 

Inclusive heavy quark production is just one example of a QCD process 
at large momentum transfer. We may consider the class of cross sections 
in which we sum over all final states that include any heavy system F = 
QQ, . . ., which can only be produced by a short-distance process in partonic 
scattering. 

Suppose for simplicity that the total mass Q of the system F is of order 
V^, the total (hadronic) center of mass energy, and that the rapidity y of 
the produced system is not a large parameter. Such a cross section can be 
expressed to leading power by the factorized expression ^ 

d'^AB-^FX 



V/ dz dXadXb(f)a/A{Xa,fJ'^)4>b/B{Xb,fJ''^ 



dQ'^dy ^ Jgys 

x6iz aab^Fx\z,e~^y—,—,Q,as{fi'^)) , (1) 

\ XaX^O j \ Xb [I j 

which is illustrated in Fig. [l|. The (/)'s are parton distributions (in some 
factorization scheme, like DIS or MS), while a is a partonic hard-scattering 
function, which at lowest order is the Born cross section for a -|- 6 -^ F -|- X, 

O = ffBorn + — <T^'^ + • • • • (2) 
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These relations may be written down for both unpolarized and polarized 
cross sections. Corrections in either case are power-suppressed 0; beginning 
at 0(1/Q^) for unpolarized or longitudinally polarized cross sections, and at 
0{1/Q) for transversely polarized cross sections. 

In Eq. (|lD, the limit Q'^/xaX^S ^ 1 is associated with corrections that 
tend to enhance the cross section. These enhancements are the object of 
threshold resummation. We will discuss below the origin of this terminology, 
why such corrections can sometimes be large, and how they can be resummed 
to all orders. Much of the following discussion follows Ref . [§ . 

The kinematics of the partonic process require that XaX^S > Q^, so that 
z < 1 in Eq. (|l]). At z = 1, the partons have just enough energy to produce 

^The extension to single-particle inclusive cross sections has recently been carried out. 



the observed final state, with no extra hadronic radiation. This is what we 
shall refer to as the elastic limit 0, or partonic threshold. It is important 
to distinguish partonic threshold from the usual concept of a threshold. In 
particular, in heavy quark production, we shall assume that the heavy quarks 
of mass Mq are produced with nonzero velocity /?, and hence with a total 
invariant mass Q^ > 4:Mq. Thus, only for (3 = does partonic threshold 
coincide with true threshold. For the Drell-Yan production of highly rela- 
tivistic lepton pairs with Q^ ^ 4m^, partonic threshold still refers to 2; = 1, 
and is a source of potentially large corrections. 
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Figure 1: Hard-scattering cross section in cut (unitarity) diagram notation. 

Typical hard- scattering cross sections a{z . . .) in Eq. (||) are distributions 
in the variable z rather than simply functions of z, because they include 
contributions from virtual as well as real gluons. We are interested in a class 
of potentially large, positive corrections due to such distributions, which 
occur in all cr*^"^ We next explain in what sense they are "large", and why 
they are positive to all orders. 

At order a", the leading logarithmic distributions in Eq. (|T|) are of the 
form (suppressing color factors, which depend on parton type) 
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whose integral with a smooth function J-'{z) (such as the convolution of 



parton distributions in Eq. (|I])) is defined as 

nl Jo 1 — z ^ 
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where we have kept only the first term in the expansion of J-'{z) about z = 
1. It is evident that such terms give, at least formally, contributions that 
grow even faster than n! at nth order. If they had alternating signs, these 
contributions might add up to a finite number somehow, but they are all of 
the same sign. 

The signs in Eq. (^) come directly from the manner in which hard-scatter- 
ing cross sections are computed. The fully inclusive Drell-Yan cross section 
da/dQ"^ illustrates the situation. The convolution in Eq. (|T|) factors into 
simple products of functions for each pair of initial-state partons a and b, 
under moments with respect to r = Q'^/S, 

^ab{N,Q) = J^ drr —^^r- 

= Y. ^a/A{N, /i) a^^{N, Q, /i) 0,/B(iV, /i) . (5) 

ab 

Identifying a'^^ as the contribution to the cross section from the parton 
combination a + b, we find 

^.7(iV,Q,/^)= n r,,, ' , ,1 tdrr^-'^^f^. (6) 

. , \ L^ rl-r tN rh- Or ii\ Jo dO/' 
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Neglecting parton labels, this is the ratio of the moment of a physical cross 
section to a product of moments of parton distributions. Because a is, by 
construction, dependent only on short-distance behavior, the ratio may be 
computed in perturbation theory, as illustrated schematically in Fig. ^ (for 
the DIS scheme). The numerator is a moment of the perturbative partonic 
Drell-Yan cross section, while the denominator is the product of moments 
of two perturbative parton distributions. For quark- ant iquark processes, the 
parton distributions are the same, so the denominator is the square of squared 
partonic amplitudes, summed over final states. 
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+ 





+ 




Figure 2: Schematic representation of moments of the Drell-Yan partonic 
hard-scattering function. 



At each order, both the numerator and denominator in Fig. ^have double- 
logarithmic terms like Eq. (^). Before moments, the perturbative Drell-Yan 
cross sections include distributions in 1 — z, like Eq. (|^) above, while in 
the deeply inelastic scattering cross section the same sort of terms depend 
on Bjorken x through 1 — x. After moments, both give double-logarithmic 
a" In^"" A^ at nth order, with A^ the moment variable. These leading loga- 
rithms are the finite remainders of corrections from n pairs of real and virtual 
gluons that attach to the scattered quarks in DIS and the annihilating pair 
in DY. In the denominator, each DIS parton distribution, which is itself of 
the form of a cross section, has both incoming and outgoing quarks, while in 
the numerator, DY involves incoming quarks only. 

On physical grounds, the numerator must vanish as 2; — > 1 and the de- 
nominator must vanish as a; — »• 1. The reason for this may be seen by recalling 
the relations of z and x to the invariant mass W of hadrons in the final state 
for the two cases: 



DY 
DIS 



W'<^Q\l-zf, 



W 



Q'^{l-x)/x. 



(7) 



The limits z ^ 1 and x -^ 1 thus both correspond to nearly elastic scattering 
0: for Drell-Yan, the annihilation of a quark pair into an electroweak vector 
boson, for DIS, the scattering of a quark into a nearly massless jet of particles. 
It is not difficult to verify that in these limits, the partonic cross section is 
suppressed |]10[- Indeed, in these limits, we expect the coherent scattering 



of hadronic bound states, whose contributions are normally suppressed by a 
power oiQ compared to incoherent partonic scattering, to dominate [11]. The 
outgoing quarks give an extra suppression in DIS when the hard-scattering 
function a, Fig. ||, is computed in perturbation theory. That is, the product 
of DIS denominators is suppressed even more than the single DY numerator. 
Then the ratio grows with moment N, from the elastic limit in z space. This 
is the source of the terms shown in Eq. (^ in a, and is the reason why they 
all have the same sign. 

The Drell-Yan cross section is the benchmark for the resummation of 
singular distributions. As above, singular distributions at z = 1 translate into 
logarithms of the moment variable N. Logarithms of A^ (to all logarithmic 
order, not just leading or next-to-leading logarithm) in the moments of the 
inclusive Drell-Yan cross section exponentiate 0, [T( 



a . {N,Q,IJ,) = CrBorn(Q) 6 



C(as)+E(N,Q,^i,as) 



(8) 



where as stands for asifi"^)- In the exponent, the function C is known to two 
loops, while the function E, which organizes all logs of A^, has the following 
form in the MS scheme, 
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The functions gi and (72 are finite series in as [0], given for MS by 
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Eqs. (|]) and (^ resum all logarithms of A^ in the sense of an order-by-order 
expansion, by reexpanding the running couplings in terms of a^. The re- 
summed integrals, however, are ill-defined for a; — > 1, no matter how large 
Q^ is, since the perturbative running coupling diverges at /x'^ = A^. Such a 
divergence is sometimes called an infrared renormalon, and may give hints 
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on the structure of power corrections to Eq. (|1]). Most importantly, it is 
possible to regulate the singularity in the running coupling without affecting 
the leading power behavior of the resummed cross section []T3||. Different ap- 



proaches to this problem can give somewhat different predictions |TJ, |T^, p!q] , 
but in any reasonable prescription the threshold corrections remain relatively 
modest, although nonnegligible in at least some physically relevant regions. 
It has been noted in several phenomenological applications that threshold 
resummation, and even fixed-order expansions based upon it, significantly 
reduces sensitivity to the factorization scale p|, |, |16|, ^, |18|, [l^ . The reason 



for this reduction is easy to see in the dependence on the factorization scale, 
/i in Eq. (|^). Indeed, for the Drell-Yan, and all related cross sections, the 
dependence on n in the hard-scattering function matches that part of the 
scale dependence of the MS distributions due to the singular part of the 
splitting functions Pii{z). That is, gi in Eq. ( [TI1| ) is exactly the coefficient of 
1/(1 — 2:)+ in Pqq[z). Because resummation matches phase space to all orders 
at partonic threshold |jlO|, all associated factorization scale dependence is also 
matched. To be specific, we may write the moments of the Drell-Yan cross 
section in a form that resums the singular threshold distributions, 

= E <t^,/A{N, /i) M,(/i, Q)crjf (iV, Q, Q) ^^/^(iV, ^^) M^ifi, Q) 

+0{l/N) , (12) 

where in the second line the factors Mg{^, Q) absorb the //-dependence in 
the exponentials of Eq. (H). These factors compensate for the singular part 
of the evolution of the parton distributions, and the //-dependence of the 
resummed expression is suppressed by a power of the moment variable. 



d 
dfi 



0,/^(iV,/i)M,(/i,Q) =Oil/N). (13) 



Of course, the importance of the remaining sensitivity to /i depends on the 
kinematics. 



3 Resummation with Color Exchange 

Beyond leading logarithms, there are important differences between the electro- 
weak-induced Drell Yan cross sections and the QCD-induced top or jet cross 
sections. These are due primarily to the presence of final-state radiation from 
scattered quarks in the latter case, which is absent in the former, and to the 
interplay of color exchange in the hard scattering with the soft radiation. As 
in the case of Eq. (^, resummation is based first of all on the factorization 
properties of the cross section in the neighborhood of the elastic limit [|I^]. 
The situation is illustrated in Fig. ^. Near the elastic limit, all gluons emitted 
into the final state have energies limited by (1 — z)Q <^ Q. Correspondingly, 
gluons with energies of order Q can appear only in virtual diagrams. Stan- 
dard factorization methods may then be used to separate the relatively soft 
but still perturbative gluons from the underlying hard scattering. This may 
be done order-by-order in perturbation theory, keeping both the hard and 
soft factors of a free of soft and coUinear divergences. The process of factor- 
ization may be thought of as the construction of an effective field theory for 
soft gluons in the presence of the hard scattering |]I|, ^ . 

In the relevant effective field theory, the incoming partons that annihilate 
into the heavy quarks and the outgoing heavy quarks themselves are rep- 
resented by ordered exponentials (Wilson lines). The Wilson lines are tied 
together in the amplitude and its complex conjugate at local vertices, Tj and 
Tj in Fig. ^, which describe the flow of color between the initial and final 
states. Indices / and J label matrices in color space. We will consider here 
the annihilation of light quarks (color indices ai and 02) into heavy quarks 
(indices 03 and 04), 

Qai (Pa) + qa2 (Pb) ^ Qa^ (Pl) + Qa, (^2) , (14) 

with kinematic invariants, 

h = {Pa ~ P2f - m'^ , ui = (pb ~ P2f ~ m'^ , s = {pa+Pbf- (15) 

In this case a convenient basis for the T's is one that represents color singlet 
and octet exchange in the s-channel. 



Other bases, particularly singlet exchange in the s- and u- channels, are also 
interesting |jl|, |^. As in Fig. ||, each choice of effective vertices leads to a 
separate soft function Sjj, which depends on (1 — z)Q only, rather than Q 
itself. 

The (virtual) hard-scattering functions, hj{Q) and hj{Q), are labelled by 
the same color exchange indices. As in most factorizations and constructions 
of effective field theories, the new vertices require renormalization. Thus we 



renormalize the soft functions M, 20 



47^ ((1 - ^)Q) = Z\j,S^;,p ((1 - z)Q) Zj,j , (17) 

and the hard functions 

h'r\Q) = ZTjh^r\Q), (18) 

where the Zkl are cutoff-dependent renormalization constants. 

Solutions to the renormahzation group equation for Su that follows from 
{ ^n\ ) are generally ordered exponentials |^, ^, but at leading logarithm in 



Sij, which is next-to- leading logarithm in the overall cross section, we can di- 
agonalize the anomalous dimension, and separate the evolution of particular 
linear combinations of composite color vertices [||. 

We can apply these results to the production of a pair of heavy quarks 
with total invariant mass Q > 2Mq = 2m at fixed rapidity y. The cross 
section is, as usual, a convolution of hard-scattering functions Cab with parton 
distributions (j)q/A and 4>q/B, as in Eq. (||), with F = QQ. Corresponding 
to the Drell-Yan result, (|^), we now have, to next-to-leading logarithmic 
accuracy, choosing fj, = Q, 

^,^QQiN) = E Sf} hj{Q)h}{Q) e^^.(-)+^- W-) , (19) 

ij 

where again as stands for as(Q^). The function C is known to one loop 
only at this time. To next-to-leading log, we need only the lowest-order 
soft functions, Sjj ~ 5jj. The function Ejj, which contains the logs of the 
moment variable A^ has a form very similar to the Drell-Yan case, but now 
with a dependence on the effective color vertices, through a third function, 
5-3, 



Jo 1-X lJ(l-x)2Q2 u'2 V 



dfi / 12 

/{l-a)2Q2 ^' 



■w/ 



= h^(Q)h*(Q) S^^ 



iJ 




Figure 3: Representation of the factorization of the hard scattering function 
a near the elastic hmit. The second part shows the soft-gluon matrix Sjj 
as a cut diagram for the scattering of incoming ordered exponentials (double 
lines - the incoming partons in the eikonal approximation) to give outgoing 
ordered exponentials (bold lines - the outgoing heavy quarks in the eikonal 
approximation). For simplicity, only a few of the possible gluon interactions 
with the ordered exponentials are shown. 



+g2{a,[Q'])+gP{a,[{l-xyQ']) 



+gi'^* («, [(1 - x)'Q 



(20) 



As before, the Qi, z = 1, 2, 3 are finite functions of their arguments. In the 
MS scheme, gi and g2 are given for incoming light quarks by ([lOD above. 
Dependence on color exchange in the hard scattering is contained entirely 
in the new functions g^ . To determine g^ , we go to a color basis that 
diagonalizes the renormalization matrix Zjj in eqs. (|I^) and (|T^. In this 
basis, 

^f H = -A,H, (21) 

where the eigenvalues A/ of the corresponding anomalous dimension matrix 
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are complex in general, and may depend on the directions of the incoming 
and outgoing partons. 

The anomalous dimension matrix of the effective vertices Tj in Fig. |^ for 
light to heavy quark annihilation in the singlet-octet basis of Eq. ([T6| ) ||l| is: 



r., 



vr 



/ -Cf{L^ + 1 + ttz) 




'^F In f"0 
Ca ^" V h ) 




\ 


21^ (it) 

\ ^ 2 


C^[41n(f^) -Lp-l-m 
-31n(f^)-ln(g)+L, 


+ Tii 


J 



(22) 

Here Lp is the vertex function in the eikonal approximation for the production 
of a pair of heavy quarks with center of mass velocity /?, 



2mVs /, 1 
- m 



/5 



P 



1 + 13 



+ ITT 



P 



Am?/ s . 



(23) 



Solving for the eigenvalues, substituting them in Eq. ([19|), and expanding the 
result to first order in a^, we can derive an explicit one- loop expression for 
the cross section for heavy quark production through light quark annihilation 
[|^, 0, ^. This result is consistent [|5[ with the explicit one-loop formulas given 
in 



21] for the single-particle inclusive cross section. 



Much the same considerations apply to electroproduction of heavy quarks. 



It is interesting to note that for the process |18 



l*{Pa) + gi{Pb) ^ QaM) + Qa,{P2) , (24) 

as for direct photon production 0, ^, |T^ there is only a single color tensor, 
coupling the produced pair to the gluon in an octet state. As a result, there 
is only one eigenvalue, which in this case is given (in the conventions of Ref. 
i) by 



a 



A. = - ( ^ 

TT 



In 



im?s 



1 + ivr 



+ {-Ca - Cp) [Lp + 1] 



(25) 



Resolved processes in photoproduction, however, involve nontrivial color mix- 



ing [^. In the limit /5 — > 0, we can make contact with quarkonia production 
in the nonrelativistic QCD (NRQCD) formalism ||23[ . 
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4 Conclusion 

The use of resuinmation as a tool for heavy-quark and other hard-scattering 
cross sections may help to put the successes and shortcomings of fixed-order 
perturbative QCD into a larger context. The reduction of factorization scale 
dependence is an encouraging observation in this direction. It will be interest- 
ing to study as well renormalization scale independence [E3]. Other directions 



for investigation include the inversion of the moments and the closely-related 
issue of power-suppressed corrections. The development of resummed polar- 
ized cross sections, for which next-to- leading order calculations are becoming 



available [E3], is sure to be a further source of insight. 
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